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Abstract 

Light scalars may be ubiquitous in nature, and their quantum fluctuations can 
produce large non-Gaussianity in the cosmic microwave background temperature 
anisotropy. The non-Gaussianity may be accompanied with a small admixture of 
isocurvature perturbations, which often have correlations with the curvature per- 
turbations. We present a general method to calculate the non-Gaussianity in the 
adiabatic and isocurvature perturbations with and without correlations, and see 
how it works in several explicit examples. We also show that they leave distinct 
signatures on the bispectrum of the cosmic microwave background temperature fluc- 
tuations. 



1 Introduction 



The inflationary paradigm has received strong observational support especially from the 
WMAP observation of the cosmic microwave background (CMB) [Ij; in a simple class of 
the single-field inflation models, density fluctuations produced by an inflaton are known 
to be nearly scale-invariant, adiabatic and Gaussian, and these properties are found to be 
in a good agreement with the observation. 

In the standard picture, the inflationary expansion is driven by a scalar field 0, the 
inflaton, which slow-rolls on a very flat scalar potential. During inflation the inflaton 
acquires quantum fluctuations, which result in slight differences in the subsequent evolu- 
tion at different places in the Universe. After inflation those differences turn into spatial 
inhomogeneities in the energy density. The nearly scale-invariant, adiabatic and Gaussian 
density perturbations are often regarded as the standard lore of the inflation theory. How- 
ever, we would like to emphasize that the lore is based on a simple but crude assumption 
that it is only the inflaton that acquires sizable quantum fluctuations during inflation. Its 
apparent success does not necessarily mean that such a non-trivial condition is commonly 
met in the landscape of the inflation theory. 

In fact, there are many flat directions in a supersymmetric (SUSY) theory and the 
string theory, and it may be even natural to expect that some of them remain light during 
inflation. If this is the case, such light scalars acquire quantum fluctuations, which may 
leave their traces in the CMB anisotropy such as isocurvature perturbations. This is 
indeed the case if those light scalars participate in the production of dark matter and/or 
baryon asymmetry of the Universe. For instance, the QCD axion |2], which was proposed 
to solve the strong CP problem, is a plausible candidate for the cold dark matter (CDM), 
and it is known that the axion generally acquires quantum fluctuations leading to the 
the CDM isocurvature perturbations. Also there are baryogenesis scenarios that contain 
light scalars. In the Affleck-Dine (AD) mechanism [3], for instance, a phase component of 
the AD field remains flat in most inflation models, leading to the baryonic isocurvature 
perturbations [U [5l [6l [7] . Although the observed density perturbation is almost adiabatic 
and no sizable isocurvature perturbation has been discovered so far, a small admixture of 
the isocurvature perturbations is still allowed. 

Non-Gaussianities have recently attracted much attention since Yadav and Wandelt 
claimed an evidence of the significant non-Gaussianity in the CMB anisotropy data [5]. 
On the other hand, the latest WMAP five-year result was shown to be consistent with the 
vanishing non-Gaussianity pTj , although the likelihood distribution of the WMAP result 
seems to favor some amount of non-Gaussianity. Also there are active studies searching 
for the non-Gaussianity [H], and it is not settled yet whether the non-Gaussianity exists 
or not. At the present stage, it is fair to say that the observations are consistent with the 
nearly scale invariant and pure adiabatic perturbations with Gaussian statistics, while 
there is a hint of non-Gaussianity at the two sigma level. 

Suppose that there is indeed significant amount of non-Gaussianity. Since it is known 
that the slow-roll inflation with a canonical kinetic term generally predicts a negligible 
amount of non-Gaussianity [IHl [HI [121 [13], we need to go beyond the simplest class of 
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inflation models. A simple and even plausible way is to introduce additional light scalars. 
In the curvaton [HI [TSl [TB] and/or ungaussiton [T7| scenarios, those hght scalars can 
generate sizable non-Gaussianity [TU [191 120] . In the presence of additional light scalars 
with quantum fluctuations, it is generically expected that isocurvature perturbation may 
arise. In particular, the non-Gaussianity hinted by the recent observations may originate 
from a small admixture of isocurvature perturbations fTEl Ul] f22\. 

We recently presented a formulation on non-Gaussianity in the isocurvature perturba- 
tions, and studied in detail how it exhibits itself in the CMB temperature anisotropy |23j . 
In Ref. [23], we found that the non-Gaussianity in the isocurvature perturbations leave 
distinctive signatures in the CMB; the non-Gaussianity is enhanced at large scales. Such 
features will enable us to distinguish the non-Gaussianity in isocurvature perturbations 
from that mainly in the adiabatic perturbation. As an example we considered the non- 
Gaussianity in the CDM isocurvature perturbations [23] and the baryonic isocurvature 
perturbations [24j. 

In this paper we extend our previous study in order to include possible correlations 
between the curvature and the isocurvature perturbations. We will give explicit examples 
in which there actually exist such correlations. Furthermore, as we did in the previous 
paper, we will present how the CMB temperature anisotropics are affected by the presence 
of the non-Gaussianities in the isocurvature perturbations correlated with the curvature 
perturbations. 

This paper is organized as follows. In Sec. [2] we extend our formalism to include 
correlation of adiabatic and isocurvature perturbations. In Sec.[3]this formalism is applied 
to some explicit models. We study features in the bispectrum of CMB anisotropy in Sec. |H 
Sec. [5] is devoted to discussion and conclusions. 

2 Formalism 

In this section, we extend the formalism developed in [23j, where the formulation to 
calculate the non-Gaussianity of the isocurvature perturbation was provided, to include 
more general case that isocurvature perturbations have correlations with adiabatic one. 
To be definite, we consider CDM isocurvature perturbation, but an extension to other 
types of the isocurvature perturbations is straightforward. 

2.1 Non-linear isocurvature perturbations and constraints 

We write the perturbed spacetime metric as 

ds^ = -U^dt^ + a^{t)e^^y,j [dx' + (5'dt) {dx^ + [^^dt) , (1) 

where TV is the lapse function, /3j the shift vector, the spatial metric, a{t) the back- 
ground scale factor, and tp the curvature perturbation. On sufficiently large spatial scales, 
the curvature perturbation if) on an arbitrary slicing at t = is expressed by [25] 

^(t/, x) = N{tf, tf, x) - log (2) 
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where the initial shcing at t = ti is chosen in such a way that the curvature perturbations 
vanish (flat slicing). Here N{tf,ti] x) is the local e-folding number, given by the integral 
of the local expansion along the worldline x = const, from t = ti to t = tf. The curvature 
perturbation evaluated on the uniform- density slicing, where the total energy density is 
spatially uniform, is denoted by (. Similarly we define the quantity d as the curvature 
perturbation evaluated on the slice where the energy density of the i-th component be- 
comes uniform {6pi{x) = 0). Hereafter we take the uniform density slicing at t = t/ in 
the last radiation dominated epoch before relevant cosmological scales enter the horizon. 
In the radiation dominated epoch the curvature perturbation is approximately given by 
( ~ (r, where Cr denotes the curvature perturbation on a slice where the energy density 
of the total radiation is spatially uniform. 

Let us assume that CDM has isocurvature fluctuation S. We allow CDM to be com- 
posed of multiple particle species each of which can have different origin. Therefore, for 
instance, some components of CDM may have isocurvature fluctuations while the remain- 
ing ones do not. 

We define the CDM isocurvature peturbation as [261 123] 

^ = 3(CcDM-Cr). (3) 

Here Ccdm is the curvature perturbation on a slice where the CDM density becomes 
spatially uniform. Then the total curvature perturbation in the matter dominated era is 
given by 

= C + (4) 

where ( is the total curvature perturbation in the radiation dominated era. 

We consider a class of models in which ( and 5* originate from the quantum fluctuations 
{S(t>a} of light scalar fields {4>a} during inflation. Note that the inflaton is also included 
in {(pa}- We can expand ( and 5* in terms of Scpa 

C = Na6(Pa + ^NabS(paS^b + • • • , (5) 
S = Sa6(f)a + ^SabS4>aS(l)b + (6) 

where A^^ = dN/d(pa and Sa = dS/dcpa B 

and summation over the repeated indices a,b, . . . 
is implicitly taken. Here we truncate the expansion at the second order and neglect higer 
order terms. For simplicity, we assume that the masses of {4>a} are negligibly small 
during and after inflation, and the fluctuations are independent to each other. Then the 
correlation functions are given by the following form, 

^ The sign of S given by Eq. ^ is opposite to the one in [25] , 

^ The evohition of {0a} is assumed to be smooth enough so that such expansions are justified. 

^ Precisely speaking, {0a} denote the field values during inflation, and they may take different values 
at the onset of oscillations. If the potentials of {0a} are significantly deviated from quadratic one, such 
differences can be important for evaluating the non-Gaussianity [19|, I29j. In this paper we neglect such 
effects. 
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with 

psm - (8) 

where k denotes the comoving wavenumber, and i^inf is the Hubble parameter during 
inflation. For later use, we also define the following: 



Before formulating the bispectra from non-Gaussianities in the isocurvature pertur- 
bations, here we define the auto/cross-correlation functions (and their spectra) of the 
primordial curvature and isocurvature perturbations, and S^. We define 

(C.^,Cg,) = (27r)='5(fci + fc2)Pc(fci), (10) 
(C,^^5,,^) = {2ixf5{k, + h)P^s{k,), (11) 

{S^,^S^,,)^{2^f5{h + h)Ps{k,). (12) 
Substituting ^ and ^ into these equations, we obtain 

1 r d^k' 

P^{k) = NaNaPs^ik) + -KbNab J j^Ps4k')Ps^{\k - k'\), (13) 

1 r d^k' 

P<:s{k) = NaSaPs^ik) + -KtSab J -^^P5^{k')Ps^{\k - k'\), (14) 

1 /■ d'^k' 

Ps{k) = SaSaPs^ + -SabSab J J^,P5^ik')Ps^{\k - k'\). (15) 

After performing the integration, the spectra P,^, P(^s and Ps can be expressed as 

P^{k) ~ [N,N, + N,,N,,^]^ ln{kL)]Ps4k), (16) 

P^sik) ~ [NaSa + NabSabA^S^P HkL)]Psm. (17) 
Psik) ^ [SaSa + SabSab^'sMf'L)]PsM- (l^) 

Here we have introduced an infrared cutoff L, which is taken to be of order of the present 
Hubble horizon scale [271 12H] • 

We define a cross-correlation coefficient by 7, 

(19) 

Uncorrelated isocurvature perturbation corresponds to 7 = and totally ( ant i-) correlated 
one is 7 = (— )1. The initial condition for the structure formation is almost adiabatic, and 
the amplitude of isocurvature perturbations is now constrained from various cosmological 
observations. For the uncorrelated case, the WMAP5 constraint is yLj 

^ < 0.190, (20) 
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while for the totaUy anti-correlated case, the constraint is 

< 0.0111. (21) 



As we will see, these constraints give upper bounds on the non-linearity parameters asso- 
ciated with the isocurvature perturbations. 

2.2 Non-Gaussianity from isocurvature perturbations 

Isocurvature perturbations must have negligible contribution to the power spectrum of the 
total curvature perturbation from observations. Therefore we can approximately obtain 

(CC) ^ (^^^^^^.) = + ^2)n(^i), (22) 

where 

P<:{k)=a<:Ps^{k), (23) 

and 

«C = + NabNabA'^s^ \n{kL). (24) 

Meanwhile, the isocurvature perturbation may significantly contribute to the three-point 
function of C^^- We define the bispectrum of C^°) B^^^ by the following equation, 

(CCC) ^ (27r)=^5(^i + k, + h) Bf^(k„ k„ k,). (25) 

This contains four kind of terms, like (CCC)) {((S), {(SS) and (SSS). For example, fo- 
cusing on the first contribution from (CCC); the bispectrum includes the following terms, 

^(A;i, k2, ks) DN^N^N^b [P54>{ki)P54>{k2) + 2 perms.] 

r d?k' - - - - (26) 

+ NabNbcN,aJ j^Ps4>{k')Ps<t>{\k, - k'\)Ps<t>{\k2 - k'\), 

and similar expressions hold for other three contributions. Performing the integration, we 
can express in terms of the four combinations of the bispectrum of C and -S", 

Bf^'iki, k2, ks) ~ (^Pccc + l^as + IPcss + ^Psss^ [Ps^{ki)Ps4.{k2) + 2 perms.] , (27) 

where each term on the right hand side (RHS) is given by 

= N^NbNab + NabNb,N,,Al^ \n{kbL) , (28) 
Pc(S = NaNbSab + 2NabNaSb + SNabNbcS^aAj^ ln( A;,L) , (29) 
Pcss = SaSbNab + 2SabSaNb + SSabSbcN^A^g^ HhL), (30) 

(^SSS = SaSbSab + SabSbcScaA'jj^ \n{kbL), (31) 
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in a squeezed configuration that one of tfie tliree wave vectors is mucfi smaller tlian tfie 
otlier two (e.g. ki <C ^2, k^) and we liave defined kb = minjfci, k2, k^}. 

In many literatures concerning the non-Gaussianity of the primordial fluctuations, 
where the adiabaticity is implicitly assumed, the magnitude of the non-Gaussianity of 
the curvature perturbation is conventionally parametrized by the so-called non-linearity 
parameter /nl which is defined by the ratio of the bispectrum to the square of the power 
spectrum. Since a single parameter is not enough to parametrize the non-Gaussianity 
in the presence of both the adiabatic and isocurvature perturbations because of their 
different effects on matter spectrum or temperature anisotropy of the CMB, we define 
four types of non-linerity parameters as following, 

6 „(cor2) _ ^-2n 6 (jso) 1 -2 n fnc)\ 

These formulae reproduce the known results for the curvaton, ungaussiton and axion, if 
only one of these non-linearity parameters exists, as we will see. However, in general, all 
the four non-linearity parameters can be concomitant and their effects on the temperature 
anisotropy have not been investigated in the previous literatures. 

3 Applications 

In this section, as an application of the formalism given in the previous section, we 
consider some simple models where the infiaton and another light field a contributes to 
the adiabatic and isocurvature perturbations. We neglect the non-Gaussianity generated 
by the infiaton itself (i.e. we set A^^^ = 0). 

3.1 Simple examples 

First we demonstrate that our formulae given in the previous section correctly reproduce 
known results for some simple models. 

3.1.1 Curvaton model with no CDM isocurvature perturvation 

If all the CDM is generated after the decay of the curvaton, no isocurvature perturbations 
are generated. Assuming No- ^ A*"^, i.e., the curvature perturbations are dominantly 
generated by the curvaton, we have 

pl^t^ = -^M^^^^ + ^L^l HhL)]. (33) 

This is the standard result for the curvaton model [I9]. But in general cases, curvaton 
models predict the existence of correlated CDM isocurvature perturbation, and it can 
significantly modify the feature of non-Gaussianity in the CMB anisotropy, as we will see 
later. 
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3.1.2 Ungaussiton model with no CDM isocurvature perturvation 

Similar to the curvaton case, but a is assumed to only affect the bispectrum of the 
curvature perturbation, and have only small effect on the power spectrum. In the limit 
A'^,^ <^ A^<^, there are no isocurvature perturbations if all the CDM is generated from the 
inflaton decay products. In this case we obtain 

^/irf = -^MN.. + NlAl HhL)l (34) 

as shown in [TTJ [20] . 

3.1.3 Axion 

Let us assume that a is stable and contributes to some fraction of the CDM. The axion 
(a) [2] is one of the well motivated candidates of such a scalar field. Since the axion has 
only negligible energy density in the radiation dominated phase, we have N„ ~ ^ 5*^ 
and (3(^(^3 ~ P^ss ~ 0, corresponding to an uncorrelated isocurvature perturbation. Thus 
we obtain 

If^T^ = ^[^'^- + Sl^l HhL)]. (35) 
This coincides with our previous result [23] . 

3.2 General curvaton model 

Let us consider a case that both radiation and CDM originate from the decay of both 
the inflaton (0) and the curvaton (o")0 The radiation and the CDM produced from the 
inflaton have the same fluctuations, i.e. the energy density of the radiation becomes 
spatially uniform on the slice where that of the CDM becomes spatially uniform. We 
denote the curvature perturbation on this slice by Ct/i- Similarly, we define as the 
curvature perturbation on a slicing on which the radiation and the CDM produced by the 
curvaton become spatially uniform. 

3.2.1 CDM from the inflaton and direct decay of the curvaton 

Let us first consider the case that some fraction of the CDM are produced directly from 
the decay of the curvaton. We assume that the universe is dominated by the radiation 
at the curvaton decays, and so we can approximate the total curvature perturbation as 
C ~ Cj.. Also we assume that the curvaton decays instantaneously when H = T^, where 
Tcr is decay rate of the curvaton. Then taking a uniform density slicing just before the 
curvaton decays, we have a relation given by 

pfix) + fpaix) = Pr{x), (36) 

The following discussions do not depend on which of them dominates the total curvature perturba- 
tion, but conveniently we call a 'curvaton' in both cases. 
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where / is the fraction of the curvaton energy density that transfers to the radiation. 
Because pf and p^- have different origin, they are inhomogeneous in general on this shce. 
Using equations pf = e^*^'''^~'')pf , pa = e^'^^''~^^pa, pr = pre^'^'^''~'^'> and C — Cr, Eq- dMD can 
be written a^ 

(1 _ e^)e4(C^-C'-) + e.e^^^^--^'-) = 1, (37) 

where = /Pc/Prldecay denotes the fraction of the radiation produced from the curvaton. 
Hence 1 — is the fraction of the radiation from the inflaton. 

In the same way, we obtain the similar equation to Eq. fl37j) for the CDM as 



(1 - ecDM)e='('^^-'^^°^) + ecDMe=^(^--^°°^) = 1, (38) 

where ecDM = (1 ~ /)Pcr/pcDM|decay denotes the fraction of the CDM produced from the 
curvaton. Eqs. (1371) and (!38l) give the fully non-linear relations between {C^r, Ccdm) and 
(C<^, Co-) under the sudden decay approximation. 

Expanding the perturbation variables in Eqs. (1371) and (!38|) up to the second order, 
we can explicitly express Cr and Ccdm in terms of C</. and Co- as 

C = (1 - R)C^ + RC. + ^i?(l - i?)(3 + R){C<t> - C.)\ (39) 

3 2 
Ccdm = (1 ~ ^cv)m)C(I> + ecDMCo- + 2^cdm(1 — ecDM)(C9i ~ Co-) , (40) 

where we have introduced a new parameter R defined by 

Meanwhile, the curvaton energy density on the uniform density slice can be written 
as Pa{x) = e^(^'^~'')po-. If we take this slice just after the curvaton starts its oscillation, 
the universe is dominated by the radiation produced from the inflatoij^ and hence C = C</>- 
Denoting the density contrast of the curvaton energy density on the uniform density 
shcing by 6^ = {pa{x) - pa)/ Pa, we find 

Co = C0 + ^log(l + 5o). (42) 



^ Precisely speaking, -when the a decays into radiation, there is a subtlety in the definition of a slicing 
on -which the energy density of a is spatially uniform. The equations in the text are valid if we interpret 
the relation between on the uniform density slicing and its background value as the definition of Co- ■ 
This is because the energy density of a on the uniform density slicing is well defined even at moment a 
decays. 

^ We assume that vacuum expectation value of the curvaton during inflation is much smaller than 
Mpi, where M^y = 2.4 x 10-'^^ GeV is the reduced Planck scale. In this case, the fraction of the curvaton 
energy density at the time when the curvaton begins its oscillations is roughly given by pa / {M'^^m'^) ~ 
a'^/M^^ ^ 1. Hence the curvaton is subdominant. Note also that it does not matter whether the inflaton 
has already decayed or not when the curvaton starts its oscillations. 
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Up to the second order in the curvaton fluctuation, the above equation can be written as 



Ca^C^ + l[S^-lsl]. (43) 
Then, using the standard formula, 

where V is the potential of the inflaton, and V^^ are its first and second derivative 
with respective to </>. Noting that 5^ — 25a/a + {5a)^ /a'^, we can express ^-nd Co- in 
terms of 50 and ^cr as 

, 2ecDM J , ecDM(l — 2ecDM) / r x2 /.^x 
+ ^5a + — {5a) . (46) 

As expected, no cross terms of 64> and 6a appear in the final expressions. From these 
results, we can immediately read the expansion coefficients as 

= - ^) > 7V.. = ^(3-4i?-2i?^), (48) 

= S^^ = 0, (49) 

1 



Sa — 2(ecDM — R) — , — 2 

(T 



eoDM(l - 2ecDM) - |(3 - 4ii - 2tf ) 



We can see that except for the trivial limit R — and ecDM = 0, isocurvature perturbation 
vanishes only when R = 1 and ccdm = 1, corresponding to the case that the curvaton 
dominates the universe before it decays and all CDM is generated by the curvaton decay 
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itself. Using this result, we can derive the following relations, 



6 (adi) 
5/NL = 

6 „{corl) 
6 p{cor2) 



6 (iso) 

^JNL 



2{i + pyR 

3 



2(l+p)2i?2 

3 



(3-4i?-2/?2), 



ecDM(l — 2ecDM) + ( -ccdm ~ R] (3 — 4i? — 2_R^) 



2(l+p)2i?3 

X 

3 



AR 2i?2 
ecDM 1 3 ^ ^ 4ecDM 



2(1 +p)2i?4 



(ecDM — Ry 



R{3 -AR- 2R^ 
R 



ecDM(l - 2ecDM) - y (3 -AR- 2R^) 



(51) 

>(52) 

(53) 
(54) 



when the fluctuation of the curvaton is dominated by the linear part, that is, 5^ ~ 25a /a. 
Here p = N^/N^ = Qa'^V'^ / {AM^^V^R^) represents the ratio of the inflaton contribution 
to the total curvaure perturbation and that of the curvaton. The limit p = corresponds 
to the standard curvaton scenario. In deriving these results, we have neglected the non- 
Gaussianity from the inflaton fluctuation because it gives /nl a value of order of the 
slow-roll parameters which are much smaller than unity [lOl [TT], [121 [13] . 

Truncating the perturbative expansion in Eqs. ( H5ll and (H6l) at the linear order, we 
find 



SI 



Nl + m 



1 



9(ecDM — R)' 

(l+p)i?2 

1 



(55) 
(56) 



Note that the cross-correlation parameter 7 is reduced to —1 in the standard curvaton 
scenario. 

Let us suppose p = for simplicity and consider the case of the large non-Gaussianity 
/nl ^10. To have the large /nl, R must be much smaller than unity. From the isocur- 
vature constraints, ecDM must be very close to R and hence ecoM is also a small quantity. 
Then the term (3c_cc is the largest and the second largest term [3qqs is suppressed by ecoM 
compared to Other terms are more suppressed by the power of ecDM- Hence the 

leading contribution to /nl from the isocurvature perturbation comes from the term (3c^c^s- 



3.2.2 CDM from the inflaton and thermal bath after the curvaton decay 

Let us next consider a case that in addition to the CDM (denoted by X) produced from 
the inflaton decay, CDM (denoted by Y) are also produced thermally during the radiation 
dominated era after the curvaton decay. We assume that all the curvaton energy density 
decays into the radiation. 
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Similar to the previous case, taking the uniform density shce just before curvaton 
decays, we obtain 



(1 - e,)e^(^^-^'-) + e,e=^(^--^'-) = 1. 



(57) 



On the other hand. Taking the uniform density shce after the total CDM is created, we 
have the following relation. 



g^g3{Cx-CcDM) _j_ g^g3{Cy "Ccdm) 



(58) 



where ex / is the fraction of X/Y in the CDM. Because CDM is assumed to be composed 
of X and y, ex + ey = 1 must be satisfied. X is produced from the decay product of the 
inflaton. Hence C,x = C<t>- Meanwhile, Y is produced from the radiation which originate 
from both the inflaton and the curvaton. Hence C,y is equal to the curvature perturbation 
on the slice where the total radiation energy density becomes homogeneous, C,y = Cr- 
Eqs. fl58p and fl57l) give the fully non-linear relation in the case of the thermally produced 
CDM. Expanding these equations up to the second order, we obtain 



C. = (1 - i?)C<A + RCa + 2^(1 - ^)(3 + Rm - C.)\ 

CcDM = C<t>~ ~ ex)(C</) ~ O 

+ ii?(l - ex) [SRex + (1 - i?)(3 + R)] {(^ - Cf. 

By doing the same procedures as in the previous subsection, we arrive at 



2Re 



X 



Srr 



cr 



[3 - AR - 2R'^ - QR{1 - ex 



(59) 
(60) 

(61) 



From this we can immediately see that the isocurvature perturbation vanishes if ex = 
0, i.e., all CDM is generated thermally after the curvaton decays, as expected. Other 
quantities (iV<^, iVo-, . . . ) are same as those obtained in the previous subsection. Thus we 
can calculate non-linearity parameters as 



,.(adi) 

/nl 



1 



2{l +pfR 
3ex 



6 „(cor2) 



2{l +pYR 



{3-AR-2R^), 

3-AR-2R^ - 2R{1 - ex) 



36? 



6 (iso) 



2{1 +p)^R 



-X 



2{i + p)m 



3 - AR - 2R' - 4R{1 - ex] 



3-AR-2R^-QR{l-ex] 



Pr 



7 



94 
1 



(62) 
(63) 
(64) 
(65) 
(66) 
(67) 
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As the simplest case, let us assume p = 0. Then, from the isocurvature constraints, 
we have a upper bound on ex as 

ex < 0.035. (68) 

We find ex must be much smaller than unity. In this case, the leading contribution to 
the non-Gaussianity from the isocurvature perturbation comes from fl^l^^^ because /nl"^^^ 
and ai'e more suppressed by the power of ex- 

Through the analyses in the previous and this subsection, we have found that in 
both cases the bispectrum of the adiabatic perturbation dominantly contribute to the 
non-Gaussianity while the leading contribution from the isocurvature perturbation comes 
from /n™^'* which are suppressed by the small parameter. But this does not necessarily 
mean the isocurvature contribution is irrelevant. Since /n™^'* represents the contribution 
from the non-Gaussian fluctuation of the correlated isocurvature perturvation between 
CDM and photon, it has more drastic effects on large-scale CMB anisotropy. Actually 
it was pointed out in [23] that an isocurvature type non-Gaussianity of the large scale 
temperature anisotropy can be about 100 times larger than the usual adiabatic case, even 
if /nl!''' = /nl°''- Thus it may be the case that the non-Gaussianity in the curvaton model 
is dominated by the isocurvature type, while satisfying the constraint on the magnitude 
of the isocurvature fluctuation from the power spectrum. We need numerical calculation 
to quantify the consequences of such mixed adiabatic and isocurvature non-Gaussianity. 



3.2.3 A 'realistic' example 

Here some particle physics motivated examples of the curvaton model are exhibited, and 
we show explicitly that it is natural to expect the existence of a correlated isocurvature 
perturbation in realistic curvaton models, if a large non-Gaussianity is generated by the 
curvaton. 

First note that R is given by 

r 

I 4M2 To 

where rricr is the mass of the curvaton, is the decay rate of the inflaton, T^- is the 
decay temperature of a, Tr is the reheating temperature after inflation defined by Tr = 
{10 / Ti'^ g^y^'^ ^yT ^Mpi with the effective relativistic degrees of freedom (7^,, T^sc is the tem- 
perature at which the curvaton begins to oscillate defined by Tqsc = (10/7r^(7^,)^/''y^?TvMpi, 
and (Tj is the initial amplitude of the curvaton, which is assumed to be much smaller than 
j^pilZl In order to generate sizable non-Gaussianity, R should be around 0.1. 

Note that in order that the curvaton can generate the observed magnitude of the density perturbation, 
R cannot be smaUer than 10^^. 



for nia > 
for m„ < 



(69) 
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In a SUSY theory, there exist many scalar fields and some of them may remain light 
during inflation. These light scalars are candidates for the curvaton [SU] ■ Good examples 
are the right-handed sneutrino [31] , flat direction in the minimal supersymmetric standard 
model [321 [7|, scalar partner of the axion (saxion) [33l [31], and so on. We assume p = 
and that the curvaton decays before the LSP freezes out T„ > mLSp/20, where mLSP 
denotes the mass of the lightest supersymmetric particle (LSP), which is assumed to be 
the lightest neutralino. In that case, the LSPs are produced in thermal bath after the 
curvaton decays, and they can be the dark matter if their annihilation cross section takes 
an appropriate value ((erf) ~ 3 x 10~^^ cm^s~^). These LSPs do not have isocurvature 
perturbation. However, in general, LSPs are also produced nonthermally by the gravitino 
decay, and gravitinos are produced at the reheating epoch after inflation [35l |36] . Thus 
energy density of the gravitino fluctuates in the same way as that of the inflaton, so 
does the LSP directly created by the gravitino decayj^ Thus those nonthermal LSPs have 
(correlated) isocurvature perturbations. This is exactly the case of Sec. 13.2.21 Here ex 
reads the fraction of the nonthermally produced LSP to the total LSP abundance. The 
gravitino number-to-entropy ratio (^3/2) is given by [351 Elj 

where rrig and m3/2 denote the mass of the gluino and gravitino, respectively. Thus 
the fraction of the nonthermally produced LSP in the total dark matter abundance is 
estimated as 

6^ ~ 5 X 10-^ f ^) ( \ . (71) 

VI TeV/ VIO^ GeV ) ^ ' 

Thus in order to satisfy the constraint fl68l) . we have upper bound on Tr as 

< 7 X 10^ GeV (l^I^ , (72) 

which is marginally consistent with an upper bound from the cosmological gravitino prob- 
lem [37] . On the other hand, Tr cannot be much smaller than this value, since otherwise R 
becomes too small, yielding too large non-Gaussianity, which conflicts with observations. 

Therefore, it seems natural to expect that some fraction of the LSP dark matter has 
different origin, which inevitably involves isocurvature fluctuations. In such a situation, 
an analysis usually done in many literatures assuming the existence of only an adiabatic 
non-Gaussianity is not valid. Instead we need to carefully study the effect of (correlated) 
isocurvature perturbations on the resulting non-Gaussian features in the CMB anisotropy. 

^ Since gravitinos decay well after the freezeout of the LSP, nonthermally produced LSPs by the 
gravitino decay remain decoupled with thermal bath. 
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4 CMB temperature anisotropy 

In this section we investigate the effects of the non-hnear isocurvature perturbation on 
the CMB anisotropy, following the notations used in [381 EH]- Since the adiabatic and 
isocurvature perturbations have quite different properties regarding their imprints on 
the CMB anisotropy, we need to correctly evaluate the bispectrum of the temperature 
anisotropy in the presence of both adiabatic and isocurvature non-Gaussianity. 

From temperature anisotropics originated from the adiabatic and isocurvature pertur- 
bations AT^'^'^^\n) and AT^'^°)(n) for a given direction n, we define agm by 



dn 



T 



T 



Transfer functions are defined by 



(adi) I 



ki 



and 



where O 



(iso) , 



it' 



(73) 



(74) 



(75) 



(adi/iso) 



[k) is the multipole moment of CMB temperature anisotropy from the 



adiabatic/isocurvature perturbation: 

^2-(adi/iso)C-X r ^3^ 



T 



^/(2£ + l)0 



(adi/iso) 



ik)P£ik-n). 



(76) 



Here P^s are the Legendre polynomials. From these equations, multipole moments can 



be divided as 



(adi) 



«£m ' + 4m where 



(adi) 



(iso) 



47rf 
47r/ 



d^k 

(2^ 

d^k 



'{k)Y;^{k)c, 



9Ti \k) ^emik) S;:. 



The anglar power spectrum of is calculated as 



Using ( ITHl) . we obtain 



(adi) 



2C, 



(cor) 



(iso) 



6f£'6„ 



^(adi) 
^£ 


2 




7T 


^(cor) 


2 




TT 


^(iso) 


2 

vr 



dk e [g^^f\k)) P^{k), 
rfA: A:2^?(?f (A:)<7?r^(A;)Pc5(A:), 
rfA;A;2(^Jf)(A;)) P5(A;). 



(77) 
(78) 

(79) 

(80) 
(81) 
(82) 
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Similarly, the angular bispectrum of is defined by 

Statistical isotropy divides the angular bispectrum into the following form, 

Here Q'^Xh^'' = J dn Yi^rm{n)Yi2m2{n)Yi.^rn:i{n) (Gaunt integral) and fe^^^j^g is the reduced 
bispectrum, on which we will focus in the following. 
Substituting fl77|) and (!78|) into (!83|) . we obtain 

g /"OO /"OO /»00 /"OO 

^^2^3 = — / dr r"^ dki kj / ciA;2 / dk^^ klji^{kir)ji>^{k2r)ji>,^{k:>,r) 
^ Jo -/o -/o -/o 

X drt^ ih) 9Th ( ^2 ) OtS^ ( ^3 ) /?CCC {ki,k2,h) 

+{fi'f/A^i)5'T?2^(^2)5'r£?(^3) + (cyclic with £'s)}%5(A;i, fca, ^s) 

+{fl'r/i'^(^i)fi'T?(^2)fi''f£3^(^3) + (cyclic with ts)}(3t;ss{ki, ^2, ^3) 

X [P5^(A;i)P5^(A;2) + (2 perms)] , (85) 

where je{x) is the spherical Bessel function. Here we have assumed that the final terms 
in Eqs. fl28l) - (!3T!) dominate. Otherwise, the above expression becomes more complicated. 

As we found in the general curvaton model (see Sec. 13.21) . requiring the isocurvature 
constraints Ps/P^ ^ 0.01 for totally anti-correlated case, the dominant contributions to 
the bispectrum come from the first two terms in Eq. ( |85l) . The first term, which is denoted 

by b^elf^e-i^ '^^^ be written as 



71' 
X 

48 
57r^ 



^Wa = ~3 / dr r'^ dki klg^^^\ki)ji^{kir)Ps^{ki) / dk2 klg^^f^ {k2)ji2{k2r)Ps4,{k2) 
' Jo Jo 

dks klgip^\k3)je.,{k3r) %f + (2 perms) 

1 POO poo 

dr r"^ I dki klgT^^\ki)ji^{kir)P(;{ki) / dk2 klg^^^ {k2)ji>2{k2r)Pc{k2) 
Jo Jo 

POO 

/ dk,klg^^^S\k,)je,ik,r) f^^i'^ + i2peTms). (86) 
Jo 
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The second term, which is denoted by ^^^j^g, can be written as 



, (corl) 



dr I dki klgi^l\ki)ji>^{kir)Ps^{ki) I dk2 k^g^Z' {^•2)ji>2{k2r)P&4>{k2) 



.2„(acii), 



dh klgf°\k3)ji^{k3r) %s + (2 perms) 



+ (cychc with €s) 



144 



dr 



dki klg!^^\ki)ji^{kir)P(;{ki) I dk2 l^ig)^^' {k2)ji2{^2r)Pc{k; 



POO 

X / dks klg^^^^Xkz) 3 i3i.hr) fi^T^^ + (2 perms 



2 (adi) , 



+ (cychc with i's). 



In deriving these equations, we have dropped logarithmic dependence on /ci, /c2, k^ of P^^i^ 

and The third and forth terms, 4i^2^3 ^^'^ ^hhh' '"^'^ written in the same 

way as 



(87) 



,{cor2) 



dr I dki klg!^^\ki)je^{kir)Ps^{ki) I dk2 k'^g)]^^^^ {k2)ji>2{k2r)P&^{k2) 



X 

Jo 
"432 



dks klgf^°\k3)ji^{k3r) P^^ss + (2 perms) 



+ (cychc with i's) 



dr I dki klg!^^^\ki)je^{kir)P(;{ki) I rf/cg /c25'T^r(^2)j£2(^2'")^c(^2) 



POO 

X / dks klg!^fj{k3)jt,.^{k3r) /jf™^^ + (2 perms 



.2„(iso), 



l,{iso) 

\e2e3 



3 
00 



+ (cychc with £'s). 



(88) 



^ I dr r^ j dki klg!j,f^{ki)je^{kir)Ps^{ki) I dk2 ki^g^T^ {k2)ji2{k2r)P&4,{k2) 



.2^(iso)/ 



X 

Jo 
1296 



dks klgf°\k3)ji^{k3r) Psss + (2 perms) 



dr r 



dki klg!j,l°^{ki)ji^{kir)P(;{ki) I d/ca k'^g}^ll'{k2)3e2{k2r)P(:{k2 



2„(iso), 



X 



POO 

/ dks klg!^fj{k3)je^{k3r) /jj™^ + (2 perms). 
Jo 



(89) 



We now move to see how non-Gaussianity from correlated isocurvature and adiabatic 
perturbations shows its signature on CMB bispectra bi-^i^i.^. To see this we plotted CMB 
bispectra h^^i^i^ in Fig. [1] and [2] using a modified version of the publicly available CMB 
code camb [50]. Here we assumed a flat scale- invariant SCDM model (fim = 1) and 
adopted a following set of cosmological parameters (fib = 0.05, ficdm = 0.95, iJo = 50), 
where fib and ficdm are the energy density of baryon and CDM, and Hq is the Hubble 
parameter in unit of km/s/Mpc. For simplicity in visualization, we have chosen following 
sets of (£1,^2) = (4, 6), (9, 11), (19, 21), (49, 51), (99, 101), (199, 201) and shown h^^^^^^ as a 
function of £3. Fig. [1] shows bispectra separately with of their dependences on initial per- 



turbations, hf^^^^, where we fixed = /^^cs = l^css = Psss = 1- 
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We can see from Fig. [1] that the bispectra with isocurvature, ^^-^^j^g; ^£i£2^3 '^^^ ^£ii2i3^ 

all 

show that their amplitudes are larger at large angular scales than that of the adiabatic bis- 
pectrum b^i^f^ia' This can be easily understood that the transfer function for isocurvature 
perturbation g'f^°\k) is large at large angular scales {^ < 10) than adiabatic one g^l^^\k) 
(See [23] for more detailed discussions). At small scales isocurvature perturbation tends 
to give small amplitude on CMB anisotropics. However, in some specfic configurations, 
bispectra arising from correlations of adiabatic and isocurvature give comparable or even 
larger amplitudes than the pure adiabatic bispectrum at relatively small scales. Especially, 
^tii2h large when two of £'s are in regime of acoustic oscillation (e.g. £1,^2 — 200) and 
the other ^ is small (£3 < 10), where g^^^"*' ^^^°\k) of their largest amplitudes are picked 
out. 

In Fig. [2] we plotted total bispectra hi^i^i^ with realistic values of f^^"* = 10, /n™^'' ~ 
—0.9, /nl'^^'' = 0.03, /nl""* = — 3 X 10""^, which, for example, can be realized by taking p = 
0,R = 0.1, ex = 0.03 in Eqs. (I63H66I) . We can see there are considerable differences from 
pure adiabatic bispectra, which we have plotted in Fig.[2]for reference. The most dominant 
contribution for deviating from pure adiabatic bispectra comes from the bispectra from 
the correlation of two adiabatic and one isocurvature inital perturbations, bf°^^(^^ since 

fl^^^^^ is at least order of magnitude larger than /nl"^ not to conflict with 

current constraints on isocurvature perturbations. Therefore we can conclude that the 
signature of non-Gaussianity from correlated isocurvature and adiabatic perturbations 
are found in CMB bispectra at least one of €s is small, where isocurvature perturbation 
gives large amplitude on CMB anisotropics. Also at intermediate angular scales {i ~ 100), 
where bispectra shows acoustic oscillation but still have sufficient power from isocurvature 
perturbation, we can see the oscillation of h^^i^^^ has different phase from that of pure 
adiabatic one h'f^l^^^, which would be a striking evidence for correlated isocurvature non- 
Gaussianity. 



5 Discussion and conclusions 

In this paper we have generalized our formalism provided in Ref. [23] for calculating 
non-Gaussianity, to include more general case where correlation between adiabatic and 
isocurvature perturbations exist, and shown that it can significantly affect the bispectrum 
of the CMB anisotropy. Actually in the curvaton scenario, the correlated isocurvature 
perturbation between CDM and radiation exists, unless all the dark matter arise after 
the curvaton decays. Although we have focused on the CDM isocurvature perturbation, 
similar mechanism can produce the baryonic isocurvature perturbation, if the baryon 
number is created before the curvaton decays. Therefore, the standard prediction for /nl 
in the curvaton/ungaussiton scenario correctly characterizes the non-Gaussian properties 
only when no isocurvature perturbations exist. In other words, there is a chance to probe 
the physics in the early Universe by using non-Gaussian isocurvature perturbations, if 
detected. 
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Figure 1: CMB bispectra h^^^^n^ (solid red line), hf^^J^ (dashed green line) hf^^J^ (dotted 
blue line) and b^i^^^i^ (dot-dashed magenta line). The thick (thin) lines correspond to 
positive (negative) values of bispectra. We plotted fe^i^j^s ^ function of £3 with fixing 
(4, £2) = (4, 6), (9, 11), (19, 21), (49, 51), (99, 101), (199, 201). Unobservable multipoles are 
shown in the shaded regions. 



The formulations provided in this paper can be applied to broad class of models. The 
right-handed sneutrino (N) may be light during inflation and acquire quantum fluctu- 
ations. The decay of generates lepton number, as well as some fraction of the total 
radiation. If the fluctuation of significantly contributes to the curvature perturbation, 
it leaves the correlation between the adiabatic and baryonic isocurvature perturbations. 
The AD field can have similar effect. The modulated reheating scenario [?T] also predicts 
mixture of adiabatic and isocurvature perturbations with correlations if the CDM/baryon 
is created before the inflaton decays. 

It is obvious that non-Gaussianity in the cosmological perturbations provides invalu- 
able information on the very early Universe. It would be interesting if the future detection 
of the non-Gaussainity tells us about the origin of CDM and/or baryon asymmetry of the 
Universe, through their small isocurvature perturbations. 

Note added 

While finalizing this manuscript, Ref. [12] appeared in the preprint server, which treats 
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4 h h 

Figure 2: CMB bispectra fe^i^afs (dashed green line) with f^^^^ = 10, /nl"^^^ = 
—0.9, f^l^'^^ = 0.03, = — 3 X lO""^. For reference, we also plotted b^e^f^es (sohd red 
line) . Same as in Fig. [H the thick (thin) lines correspond to positive (negative) values of 
bispectra and unobservable multipoles are shown in the shaded regions. 



similar subjects to ours. 
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